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Abstract 

We show that, contrary to recent claims, the J/'i/' suppression resulting from 
its interaction with comovers does not saturate at large transverse energy Ej-- On 
the contrary, it shows a characteristic structure - change of curvature near the knee 
of the Et distribution - which is due to the Et (or multiplicity) fluctuation, and 
agrees with recent experimental results. 
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An interesting result of the last run (1998 data) of the NA50 collaboration at 
CERN on the transverse energy (Et) dependence of J/ip suppression in PbPb collisions, 
is the observation [|[| of a convexity at large Et- More precisely, for E^ > 100 GeV 



ratio R{Et) of J/ip over Drell-Yan (DY) cross-sections increases with increasing Et- In 
sharp contrast with this result, models of the J /ip suppression in non-quark-gluon plasma 
(QGP) scenarios [2-9] - such as the one based on the interaction of the J/ip with comovers 
- exhibit a clear saturation at large Et- 

In this work, we show that the above feature of the comovers model is only true 
up to the knee of the Et distribution {Et ~ 100 GeV). Beyond this value, we enter 
into the tail of the Et distribution - where the increase in Et is due to fluctuations. 
This fluctuation, which has not been taken into account in most calculations, produces 
a corresponding increase in the density of comovers - which, in turn, increases the J/%1) 
suppression at large Et- 

In order to illustrate this phenomenon we use the model introduced in ref. 0. 
Here, as in most non-QGP models, the J/ if) suppression is due to two mechanisms : 
absorption of the pre-resonant cc pair with nucleons (the so-called nuclear absorption) and 
the interaction of the J /if) with comovers. The corresponding J/ip survival probabilities 
are given by [0]. 



(which corresponds to the so-called knee of the Et distribution [Q]), the slope of the 
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The survival probability S'^° depends on the density of comovers Ny°[b, s) in the rapidity 



region of the dimuon trigger 2.9 < yiab < 3.9 and Nf = 1.15 fm~^ |0, Q is the corresponding 
density in pp colhsions. In order to compute iV™, various hadronic models have been 
used in the hterature. For instance in ref. ||^ it has been assumed that the hadronic 
multiphcity is proportional to the number of participant nucleons (the so-called wounded 
nucleon model), while in ref. a formula based on the dual parton model (DPM) was 
used - which includes an extra term proportional to the number of binary interactions. In 
this paper we use the DPM formula (eq. (6) of [0]). In both cases, the calculations do not 
include the fluctuations mentioned above and, therefore, cannot be applied beyond the 
knee of the Ex distribution - where the increase in Ex (or multiplicity) is entirely due to 
fluctuations. In order to introduce these fluctuations, it is convenient to recall the other 
formulae needed to calculate the J/ip suppression. 

At fixed impact parameter b, the J/tp cross-section is given by 

a^^(6) = 1m I dh m{b, s) S'^'^ib, s) ^="(6, s) , (3) 

where m{b, s) = AB Cpp Ta{s) Tb(6 — s). The corresponding one for DY pair production 
is obtained from (H) putting aabs = o"co = (i.e. S'^'"^ = S*™ = 1) and is proportional 
to AB. In this way we can compute the ratio of J/ip over DY as a function of the 
impact parameter. Experimentally, however, the ratio R{Et), is given as a function of the 
transverse energy Et measured by a calorimeter, in the rapidity interval 1.1 < yiab < 2.3. 
In order to compute R{Ej) we have to know the correlation P{Et, b) between E^ and 
impact parameter, which is given by |3] 



P{Et, b) = , = exp 
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Here 



E^^ib) = q NZip) + k[A - mAmE^n , (5) 

mA{b) is the number of participants of A (at fixed impact parameter), Ein = 158 GeV/c is 
the beam energy and k = 1/4000 0. In @ and (H) N^°i{h) is obtained by integrating the 
comover density Ny°{b,s) over (Ps, and dy (in the rapidity range of the Et calorimeter). 
The second term in (^) was introduced in ref. in order to reproduce the correlation 
between Et and the energy Ezdc of the zero degree calorimeter. It was interpreted as due 
to intra-nuclear cascade - which is present here due to the location in rapidity of the NA50 
calorimeter. This term is sizable for peripheral collisions, when many spectator nucleons 
are present, and dies away for central ones. The parameters q = 0.56 and a = 0.94 are 
obtained from a fit to the minimum bias E^ distribution at large Et- The parameter q 
gives the relation between multiplicity of comovers (positive, negative and neutrals) and 
the Et of the NA50 calorimeter (which contains only neutrals). The product qa controls 
the width of the Et distribution at fixed 6 ~ 0. The J/ip and DY cross-section at fixed 
Et are then given by 



j d% afl^'^^ P{EtM ■ (6) 
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The quantity E^^{b) in eq. (^) does not contain fluctuations - hence the index 
NF. This is obvious from the fact that the parameter a is not present in (^). In order to 
see it in a more explicit way, we plot in Fig. 1 the quantity 

F{Et) = Et/E^^{Et) , (7) 



where 



We see that coincides with Et only up to the knee of the Et distribution. Beyond 
it, E^^ is smaller than the true value of E^- This difference is precisely due to the E^ 
fluctuation. 

As discussed above, in order to compute the ratio R{Et) beyond the knee of the 
Et distribution it is necessary to introduce in N^" the E^ (or multiplicity) fluctuations 
responsible for the tail of the distribution. In order to do so, we use the experimental 
observation that multiplity and Et distributions have approximately the same shape. 
This indicates that the fluctuations in Et are mainly due to fluctuations in multiplicity - 
rather than in pt- This leads to the following replacement in eq. (|^): 



iV-(6,s)-.iV-(6,s) F{Et) . (9) 

In this way the results for the ratio R{Et) are unchanged below the knee of the distribution 
(see Fig. 1). Beyond it, the J/ip suppression is increased as a result of the fluctuation. 

We turn next to the numerical results. In ref. [0] we used for the two parameters 
of the model aabs = 6.7 mb and dco = 0.6 mb. In this case, the computed J/ip suppression 
at Et ~ 100 GeV is somewhat too small 0. Clearly, we can increase it by increasing 
the value of o"co- However, we then increase the value of the suppression for peripheral 
collisions. This, in turn, can lead to some conflict with the SU data (see for a discussion 



on this point). However, recent data [jTOl on the J/ip cross-section in pA collisions, point 
to a smaller value of aabs - of 4 to 5 mb. With this value of aabs, we can increase aco from 
0.6 mb up to 1.0 mb without decreasing the J/ip suppression for peripheral collisions. In 
Fig. 2 we present the result of our calculation using aabs = 4.5 mb and aco = 1 mb. We see 



that the main features of the data are reproduced. In particular our curve shows a shght 
change of curvature at Et ~ 100 GeV, which is entirely due to the effect of fluctuations - 
and is seen in the 1998 NA50 data [|l|. The physical origin of this change in the slope of 
Ej- is the following: when approaching the knee of the E^ distribution from below, the 
number of participants approaches 2A and changes slowly. The latter is also true for the 
multiplicity of comovers. Beyond the knee, the multiplicity increases faster due to the 
fluctuations and produces a faster decrease of R{Et) (see Fig. 1). 

We want to stress that the shape of our curve in the lower half of the Et region 
(where the ratio R{Et) changes rather fast with Et) is sensitive to the relation between 
Et and impact parameter. We see from eq. (|^) that this relation depends on the size 
of the contribution of the intra- nuclear cascade (parameter k). As mentioned above, the 
value k = 1/4000 used here was obtained in [0 from the best fit of the correlation between 
Et and the energy Ezdc of the zero degree calorimeter. However, since we do not have 
a totally reliable expression for the latter, there is an uncertainty in the value of k. In 
order to illustrate its effect on R{Et), we show in Fig. 2 (dashed line) the result with 
k = 1/2000, i.e. doubling the (comparatively small) contribution of the intra-nuclear 
cascade. The effect is concentrated in the lower half of the Et range. This uncertainty 
would not be present if the J/ip suppression were given as a function of either Et or 
charged multiplicity at mid-rapidities. 

The large Et structure seen by the NA50 collaboration can also be explained as- 



suming a deconfining phase transition |[Tl[. At the energies of the Relativistic Heavy Ion 
Collider (RHIC) at Brookhaven, it will be possible to determine which of these two mech- 
anisms is the correct one. Indeed, the transverse energy (or the corresponding energy 
density) where the structure has been seen by NA50, will be reached at RHIC well below 



the knee of the Et distribution, and, if our interpretation is correct, no structure will be 
present. It will, however, appear at higher values of Et - when reaching the knee of the 
Et distribution at = 200 GeV. 
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FIGURE CAPTIONS 



FIG. 1. The ratio F{Et) in eqs. (0), (|). 



FIG. 2. The ratio R{Ej) of J/ip over DY cross-sections, obtained with aabs = 4.5 mb 



and a CO = 1 nib, compared to the NA50 data [|l| |[T^. The full curve corresponds to 
k = 1/4000 0] in eq. (|). The dashed curve is obtained with k = 1/2000 (see main text). 
The black points correspond to 1996 Pb-Pb data, the black squares correspond to 1998 
Pb-Pb data, the white points to 1996 analysis with minimum bias and the white squares 
to 1998 analysis with minimum bias. 
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